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ABSTRACT 


!.  V 

A new  image  restoration  system  that  is  applicable  to  the  problem  of 
restoring  an  image  degraded  by  blurring  and  additive  noise  is  presented, 
in  thic  report.  The  system  is  developed  by  attempting  to  estimate  more 
accurately  the  frequency  response  of  typical  image  restoration  filters 
available  in  the  literature.  The  resulting  system  combined  with  its 
short  space  implementation  computationally  simple  and  appears  to  compare 
quite  well 
are  given 


in  performance  with  other  restoration  techniques.  Some  examples 
to  illustrate  the  performance  of  the  new  image  restoration  system 


I. 


INTRODUCTION 


A number  of  techniques  (1,2)  have  been  proposed  in  the  literature  to 
restore  an  image  from  degradations  due  to  blurring  and  additive  noise. 

Wiener  filtering  (3,4),  parametric  Wiener  filtering  (2),  power  spectrum 
filtering  (5)  and  geometrical  mean  filtering  (5)  are  some  of  the  examples. 

In  many  cases,  the  restoration  techniques  are  based  on  the  assumption  that 
an  image  can  be  modelled  by  a stationary  random  field,  and  restoration  is 
achieved  by  filtering  the  degraded  image  with  a linear  space  invariant 
restoration  filter,  the  frequency  response  of  which  is  a function  of  the 
ideal  image  power  spectral  density. 

There  are  at  least  two  difficulties  in  this  approach.  For  a typical 
image,  each  part  of  an  image  generally  differs  sufficiently  from  other 
parts  so  that  the  stationarity  assumption  over  the  entire  image  is  not 
generally  valid.  The  second  difficulty  is  that  in  practice*  the  image 
power  spectral  density  is  not  given  and  has  to  be  estimated.  A common 
procedure  used  is  to  estimate  the  power  spectral  density  by  an  average 
over  different  images  or  noise-free,  blur-free  prototype  images  whose 
contents  are  similar  to  the  suspected  content  of  the  ideal  image.  Such 
a procedure  is  typically  quite  cumbersome,  and  the  resulting  power  spectral 
density  estimate  is  at  best  a crude  approximation  to  the  true  image  power 
spectral  density.  In  this  paper,  we  develop  a new  image  restoration 
technique  by  attempting  to  reduce  the  effects  of  these  problems. 

To  suggest  a basis  for  the  new  restoration  technique,  consider  an 
image  f(n^,n2)  degraded  by  additive  noise  and  suppose  that  our  approach 
is  to  restore  the  degraded  image  by  a restoration  filter  such  as  a Wiener 
filter.  To  implement  a Wiener  filter,  we  need  a good  estimate  of  P^foj^.u^) , 
the  power  spectral  density  of  f(n1,n2).  A reasonable  approach  to  obtain 
Pf(wl»u2)  is  to  estimate  it  from  f(n^,n2).  To  get  a good  estimate  of 
f(n1,n2),  however,  we  need  a good  estimate  of  Pf (coj ,ou2) . This  then  suggests 
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an  iterative  procedure  in  which  we  begin  with  an  initial  estimate  of 
P^.(oj i.tu 2)  and  then  iteratively  estimate  f(nj,n2)  and  until  a 

converging  solution  or  a desired  performance  is  achieved.  At  a first  glance, 
this  iterative  procedure  to  estimate  f(n^,n2)  does  not  appear  to  be  very 
useful,  since  many  iterations  may  be  required  to  obtain  a converging  solution 
and  the  resulting  solution  may  be  very  sensitive  to  the  initial  estimate  of 
Pf  (u)-i, tu • As  will  be  seen  later,  however,  with  a particular  method  of 
estimating  Pf(o)1,u)2)  from  the  estimated  f(n1,n2)  and  with  a particular 
choice  among  existing  restoration  filters,  a converging  solution  that  does 
not  depend  on  the  initial  estimate  of  PjCuij.u^)  can  be  obtained  analytically. 
An  image  restoration  technique  based  on  this  converging  solution  for 
f(nj,n2)  will  be  seen  to  be  computationally  simple  and  to  perform  well  as 
an  image  restoration  system. 

To  reduce  the  effect  of  the  non-stationarity  problem,  the  restoration 
technique  developed  in  this  paper  will  be  implemented  on  a short  space 
basis  in  which  an  image  is  divided  into  many  subimages  and  each  subimage 
is  restored  separately  and  then  combined.  As  will  be  discussed  later  and 
consistent  with  other  results  (1,6)  in  the  literature,  short  space  imple- 
mentation when  performed  properly  achieves  a better  system  performance  than 
the  approach  to  process  an  entire  image  at  one  time. 

The  overall  objectives  of  this  paper  are  to  develop  a new  image 
restoration  system  and  illustrate  its  performance.  In  Section  II,  we  develop 
a new  image  restoration  system.  In  Section  III,  we  discuss  various  issues 
related  to  short  space  implementation  of  the  system  developed  in  Section  II. 
In  Section  IV,  we  discuss  and  illustrate  various  examples  in  which  the  new 
restoration  system  is  applied  to  restore  images  degraded  by  additive  noise, 
or  blurring  and  additive  noise. 
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II.  IMAGE  RESTORATION  BY  SPECTRAL  SUBTRACTION  (SSIR) 


The  problem  that  we  consider  in  this  paper  is  to  restore  an 
image  degraded  by  blurring  and  additive  background  noise  as  shown  in 
Figure  1.  The  degradation  model  in  Figure  1 can  also  be  represented  by 


g(nj,n2)  = f (n1>n2)  * bfn^np  (1-a) 

y(n1,n2)  = gfoj.n^  + dfn^n^  (1-b) 

and  thus  yfn^np  = ffn^n^  * bfn^n^  + dfn^n^  , (2) 

where  f(n1>n2)  represents  a digital  image, 

b(nj,n2)  represents  a linear  space  invariant  point  spread  function, 

and  d(n1>n2)  represents  an  additive  noise  component. 

The  restoration  problem  to  be  considered  is  to  restore  f(n^,n2)  from  the 
degraded  image  y(n1>n2). 

The  restoration  problem  discussed  above  has  been  studied  extensively 

in  the  literature  and  a number  of  techniques  have  been  proposed.  Many  of 

them  can  be  viewed  as  a two  step  procedure  in  which  gfn^n^  is  estimated 

from  y(n1,n2)  by  a noise  reduction  system  and  then  f(nj,n2)  is  estimated 

from  g(nj,n2)  by  a deblurring  system.  More  specifically,  consider  the 

restoration  system  shown  in  Figure  2.  In  the  figure,  P (ioj.g^),  Bfo^.u^) 

„ 1 8 
and  g(n1,n2)  represent  the  power  spectrum  of  g(n1>n2),  the  discrete  space 

Fourier  transform1  of  b(n1>n2)  and  the  estimate  of  gO^.n^  respectively. 
Pd(ai1,w2)  and  f(nJ,n2)  are  similarly  defined,  and  a and  6 are  constants. 

If  B is  unity,  the  restoration  system  in  Figure  2 corresponds  to  the  para- 
metric Wiener  filtering  technique  (2)  and  if  a is  also  unity,  it  reduces  to 
the  Wiener  filtering  technique  (3).  If  a is  unity,  6 is  S and  BfWj.o^)  is 
zero  phase,  the  system  corresponds  to  the  power  spectrum  filtering 

1.  The  definitions  of  these  terms  are  given  in  the  Appendix.  See  also  (7). 
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technique  (5).  If  a and  8 are  some  constants,  then  the  system 
corresponds  to  the  geometrical  mean  filtering  technique  (2,5). 

In  implementing  the  noise  reduction  system  Hj(<dj,oi)2)  in  Figure  2, 

PgCajj.u^)  and  have  to  be  known^.  In  many  cases  (2),  the  estimation 

of  P^ (oij,^)  is  relatively  simple.  To  estimate  P^(ujj ,w2) , a common  pro- 
cedure used  is  to  estimate  first  by  averaging  the  spectral  density 

over  many  different  images  or  over  noise-free,  blur-free  prototype  images 
and  then  using  the  relationship  of 

pg(a>l.w2)  = | B (coj , u>2)  | 2 3 • Pf(co1,u)2)  (3) 

Since  different  images  or  even  subimages  generally  have  different  power 
spectral  density,  a more  reasonable  approach  would  be  to  estimate  P^(gjj,u>2) 
from  gOij.n,,).  The  function  gCn^.n,)  is  not  known,  of  course,  but  could  be  estimated 
by  the  noise  reduction  system  in  Figure  2 which  requires  the  knowledge  of 
Pg(u)j,u)2)*  This  then  suggests  an  iterative  procedure  in  which  we  begin  with 
an  initial  .imate  of  P^fiOj,^)  and  then  iteratively  estimate  gCn^.n^  and 
P (u)j,u»2)  until  a converging  solution  or  a desirable  performance  is  achieved. 

Since  many  iterations  may  be  required  to  obtain  a converging  solution, 
the  iterative  procedure  discussed  above  to  estimate  g(n^,n2)  appears  to  be 
computationally  undesirable.  However,  with  a particular  method  of  estimating 
Pg(ojj,(i)2)  from  g(nj,n2)  and  with  a particular  choice  of  8,  a simple  converging 
solution  can  be  analytically  obtained.  More  specifically,  suppose  we  let  the 

■I  ^ 

value  of  8 be  *5  and  estimate  P (to,  .uu)  by  r-| G (on  , uu ) I where  k is  a scaling 
factor  that  normalizes  the  power  and  energy  spectral  densities  . Then, 
from  Figure  2,  §(nj,n2)  at  the  ith  and  i+lth  iteration  is  related  to  each 
other  by 


Gi+l(u)l,u)2)  = Y (Uj  ,0)2) 


|Gi(u)1,U)2)  | 

Gi(u1,w2)|2  + a • k-  Pd(di1,w2) 


2.  In  the  literature,  there  exist  techniques  (8)  in  which  an  explicit  knowledge 
of  P (copu^)  is  not  necessary.  Even  in  such  cases,  P^fojpU^)  still  needs 

to  b§  known. 

3.  See  Appendix. 
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where  Y(u>  ,oj  ) represents  the  discrete  space  Fourier  transform  of  y(n..,n  ), 

A 1 2 A 1 dZ 

and  G.  and  G^+1  ,gj2)  represent  the  discrete  space  Fourier  transform 

of  gfn^.n^)  at  the  ith  and  i+lth  iteration.  If 


Gi+1(u)1,ai2)  = Gi(aj1,u2) 


(5) 


then  a converging  solution  of  the  iterative  algorithm  are  reached  since 
further  iterations  will  produce  identical  results.  Combining  equations  (4) 
and  (5),  a converging  solution  for  G(wi,u>2)  is  given  by 

$G(o)i,u)2)  = < Y(u)i,U)2) 

and  iGfWpU^I  = (6) 

^|y(u1,w2)  |2  - a • k • Pd(iu1,a;2)^'5  if  | Yfuij  ,u>2)  | 2 > a • k • P^ajj.u^) 

0 otherwise, 

where*  G(a)j  ,0^)  and  *¥(0^,0^)  represent  the  phase  of  g(n.j,n.,)  and  y(nj,n2). 
Equation  (6)  completely  characterizes  the  new  noise  reduction  system.  The 
phase  of  g(nj,n2)  is  estimated  by  the  phase  of  y(n1>n2)  and  the  transform 
magnitude  of  g(n^,n2)  is  estimated  by  a particular  form  of  spectral  sub- 
traction. 

Combining  the  noise  reduction  system  given  by  equation  (6)  with  the 
deblurring  system  in  Figure  2,  the  resulting  overall  image  restoration  system 
is  shown  in  Figure  3.  We'll  refer  to  this  algorithm  as  "Spectral  Subtraction 
Image  Restoration"  (SSIR)  technique.  As  is  clear  from  Figure  3,  the  SSIR 
technique  does  not  require  an  explicit  estimation  of  P^Cu^.u^)  or 
P^fu^jU^)  and  the  computations  required  for  restoring  an  image  are  not  any 
greater  than  typical  restoration  techniques  such  as  Wiener  filtering. 

The  SSIR  technique  was  developed  based  on  a converging  solution  of  an 
iterative  procedure  that  capitalizes  on  a particular  restoration  filter 
(6=1/2)  and  a particular  method  of  estimating  the  power  spectrum  01  a given 
data.  In  the  iterative  procedure,  we  could,  of  course,  have  used  a different 
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value  of  3 or  a different  method  of  estimating  the  power  spectrum  of  a given 
data.  In  such  a case,  however,  a converging  solution  cannot  generally  be 
obtained  analytically  and  the  resulting  algorithm  would  be  unattractive  from 
the  computational  point  of  view. 


III.  SHORT  SPACE  IMPLEMENTATION 

The  noise  reduction  system  developed  in  the  previous  section  is  based  on 
the  assumption  that  y(nj,n2)  can  be  approximately  modelled  by  a stationary 
random  field.  For  a typical  image,  each  part  of  an  image  generally  differs 
sufficiently  from  other  parts  so  that  the  stationarity  assumption  is  not 
generally  valid.  To  reduce  the  effect  of  the  non-stationarity  problem,  a 
reasonable  approach  is  to  implement  the  noise  reduction  system  on  a short 
space  basis  in  which  the  degraded  image  is  divided  into  many  subimages  and 
each  subimage  is  restored  separately  and  then  combined.  More  specifically, 
we  first  apply  a short  space  window  function  w^Cn^.n^  to  yf^.n^  in 
equation  (1-b)  so  that 

y(nrn2)  • = g(nj,n2)  • w^fn^n^  + d(nj,n2)  • ^ (n^)  . (7) 

Rewriting  equation  (7), 

XijCW  = «ijCni»n2)  + dijfni'n2)  (8) 

The  noise  reduction  system  developed  in  Section  II  is  then  applied  to 
equation  (8)  to  recover  g^Cn^n^  from  y^fn^np.  fKn^n^  is  then 
constructed  by  combining  ^.(n^^)  as  follows; 

g(n1,n2)  = I E (9) 

i j J 

Then  the  deblurring  system  is  applied  to  g Cn^ ,n2)  to  restore  f(nj,n2). 

To  successfully  implement  the  image  restoration  system  on  a short  space 
basis,  the  window  function  w^Cn^n^  must  be  carefully  chosen.  For  example, 
to  use  equation  (9)  in  constructing  an  image  from  its  subimages, 
w^(nj,n2)  has  to  satisfy  the  following  equation; 


E Ew^.(nj,n2)  = 1 for  all  nj,n2  of  interest 

i j ij 
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In  addition,  w^^fn^.n^)  is  desired  to  be  a smooth  function  to  avoid  some 
possible  discontinuities  that  may  appear  at  the  subimage  boundaries  in  the 
processed  image. 

One  way  to  find  a smooth  2-D  window  function  that  satisfies  equation  (10) 
is  to  form  a separable  2-D  window  from  1-D  windows  that  satisfy  similar 
conditions.  For  example,  let  wxCnj)  and  wy(n2)  represent  smooth  1-D  windows 
with  durations  of  2K  and  2L  respectively; 

w^O^)  = 0 except  for  0 <_  n^  <_  2K-1 

w (n.)  = 0 except  for  0 _<  n,  <_2L-1 

y z z 

Now  let  w^(np  and  w^  (n^)  be  defined  by 


wi(n1)  = w^n^iK)  (11-a) 

and  Wj(n2)  = wy(n2-jL).  (11-b) 

If  w^(nj)  and  w^  (n2)  satisfy 

E w.(ni)  = 1 (12-a) 

and  E w. (n  ) = 1 , (12-b) 

1 


it  is  straightforward  to  show  that  a separable  2-D  window  w^  fn^np  = 
Wjfaj)  • Wj  (n2)  is  smooth  and  also  satisfies  equation  (10).  Two  such 
window  functions  are  2-D  separable  triangular  or  Hamming  windows  over- 
lapped with  its  neighboring  window  by  half  the  window  duration  in  each 
dimension.  The  case  of  a 2-D  separable  triangular  window  is  shown  in 
Figure  4. 
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In  summary,  in  estimating  gfrj.n^  from  yO^.n^  by  the  SSIR  technique, 
we  first  divide  yfrj.n^  into  subimages  by  applying  a smooth  window  and  the 
noise  reduction  system  is  applied  separately  to  each  subimage.  The  resulting 
subimages  are  then  combined  to  form  g(nlfn2)  which  is  then  deblurred  to 
restore  f(nj,n2).  In  the  next  section,  we  illustrate  the  performance  of  the 
SSIR  technique  when  applied  to  restore  images  degraded  by  additive  noise, 
or  blurring  and  additive  noise. 


V. 


EXAMPLES  AND  DISCUSSIONS 


In  this  section,  we  illustrate  and  discuss  some  examples  where  the  short 
space  SSIR  technique  developed  in  Sections  II  and  III  is  applied  to  restore 
images  degraded  by  additive  noise,  or  degraded  by  blurring  and  additive  noise. 
We  first  consider  the  case  of  images  degraded  by  additive  noise  alone.  In 
Figure  5 are  shown  two  nois^-free,  blur-free  images  of  256x256  pixels  with 
each  pixel  represented  by  8 bits.  In  Figure  6 are  shown  two  degraded  images 
at  S/N  ratio  of  lOdb  which  were  generated  by  adding  zero-mean  white  Gaussian 
noise  to  the  images  shown  in  Figure  5.  The  Gaussian  noise  was  digitally 
generated  and  S/N  ratio  is  defined  by 

Variance  of  f(n.,n_) 

S/N  ratio  = 10  • log  

Variance  of  dfaj.nj)  (13) 

In  Figure  7 are  shown  the  images  restored  by  applying  the  short  space  SSIR 
technique  to  the  noisy  images  in  Figure  6 with  a = 2.5  and  the  subimage  size 
(window  size)  of  32  x 32  pixels. 

The  value  of  a in  the  system  was  chosen  as  a compromise  between  noise 
reduction  and  signal  distortion.  If  a is  very  small,  the  short  space  SSIR 
technique  approaches  an  identity  system.  If  a is  very  large,  then  large 
noise  reduction  can  be  achieved  but  at  the  expense  of  signal  distortion. 

Based  on  a few  other  examples  in  addition  to  what  are  shown  in  this  paper, 
we  have  found  that  a reasonable  choice  of  a at  which  substantial  noise 
reduction  is  achieved  without  noticeable  signal  distortion  is  in  the  range 
of  1.5  to  3.0. 

The  subimage  size  in  the  short  space  implementation  was  chosen  based 
on  two  considerations.  If  the  subimage  size  is  too  large,  then  the  non- 
stationarity  problem  that  we  discussed  in  Section  III  could  affect  the 
system  performance.  On  the  other  hand,  if  the  subimage  size  is  too  small, 
the  ability  of  the  system  in  reducing  noise  may  be  diminished  due  to  lack  of 
available  data.  For  the  types  of  images  that  we  considered  such  as  pictures 
of  "clock"  and  "aerial  view  of  a village",  the  subimage  size  in  the  range  of 


16  x 16  pixels  to  32  x 32  pixels  produced  the  best  results.  Consistent  with 
the  above  discussions,  when  the  SSIR  technique  was  applied  to  the  degraded 
images  in  Figure  6 with  a too  large  subimage  size  or  a too  small  subimage 
size,  the  system  performance  was  not  as  good. 

In  generating  the  images  shown  in  Figure  6,  we  used  a 2-D  separable 
triangular  window  discussed  in  Section  III.  Two  other  windows,  namely  2-D 
separable  Hamming  and  rectangular  windows  which  both  satisfy  equation  (10) 
were  also  considered.  The  results  obtained  using  a Hamming  window  were  very 
similar  to  those  of  a triangular  window.  However,  the  results  obtained  using 
a rectangular  window  had  the  discontinuity  problem  at  the  subimage  boundaries. 

Now  we  consider  the  case  when  images  are  degraded  by  both  blurring  and 
additive  noise.  In  Figure  8 are  shown  two  images  obtained  by  first  blurring 
the  images  of  Figure  5 with  a Gaussian  shaped  point  spread  function  and  then 
adding  white  Gaussian  noise  at  S/N  ratio  of  20db.  In  Figure  9 are  shown 
the  restored  images  by  the  short  space  SSIR  technique  with  a = 2.5  and  a 
separable  2-D  triangular  window. 

The  value  of  a,  subimage  size  and  window  shape  were  again  chosen  based 
on  our  previous  discussions.  In  implementing  the  deblurring  system  in  Figure  3, 
| B(gl)j  ,0)^)  I)  the  magnitude  of  Bfajj,^),  was  first  compared  to  some  threshold 
value.  If  the  magnitude  is  smaller  than  the  threshold  value,  | BCu^  | was 
set  to  the  threshold  value  without  affecting  the  phase  of  B(a)j .a^) • This  step 
was  taken  to  avoid  possible  singularity  problems  of  inverse  filtering  when 
Bf'ii^.to^)  is  close  to  zero.  In  general,  inverse  filtering  accentuates  the  back- 
ground noise  and  is  useful  only  when  the  background  noise  level  is  very  low. 
Thus,  for  example,  when  we  performed  inverse  filtering  on  the  images  in  Figure  8 
without  applying  the  noise  reduction  system,  the  resulting  images  were  dominated 
by  background  noise.  When  accompanied  by  an  effective  noise  reduction  system, 
however,  inverse  filtering  can  be  used  as  a deblurring  system  as  evidenced  in 
Figure  9. 

A serious  attempt  has  not  yet  been  made  to  compare  the  performance  of 
the  short  space  SSIR  technique  to  other  restoration  techniques  in  the 
literature.  However,  a very  informal  comparison  was  made  based  on  several 
images  recovered  by  the  restoration  techniques  that  can  be  represented  by 
Figure  2.  Preliminary  results  indicate  that  the  performance  of  the  restoration 
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Fig.  6.  Images  of  Fig.  5 degraded  by  additive  noise  at  S/N  ratio  of  10  dB. 


Fig.  8.  Images  of  Fig.  5 degraded  by  blurring  and  additive  noise  at  S/N 
ratio  of  20  dB. 


system  developed  in  this  paper  compares  quite  well  with  other  restoration 
techniques.  A more  detailed  comparison  requires  further  study. 

Finally,  in  this  paper  we  have  considered  short  space  processing  for 
image  restoration.  Even  though  the  notion  of  short  space  processing  was 
considered  in  the  context  of  a specific  image  restoration  technique,  it  can 
also  be  applied  to  other  restoration  techniques.  Some  preliminary  results 
indicate  that  short  space  processing  based  on  smooth  2-D  windows  such  as  a 
separable  triangular  or  Hamming  window  can  improve  the  performance  of  other 
restoration  techniques  such  as  Wiener  filtering. 

V.  CONCLUSION 

In  this  paper,  we  developed  a new  image  restoration  system  which  is 
applicable  to  the  problem  of  restoring  an  image  degraded  by  blurring  and 
additive  noise.  The  restoration  system  consists  of  a noise  reduction  system 
and  a deblurring  system.  The  main  features  of  the  noise  reduction  system  are 
to  estimate  the  spectral  magnitude  of  the  blurred  image  by  a particular  form 
of  spectral  subtraction  and  to  implement  the  system  on  a short  space  basis. 
The  deblurring  system  is  a modification  of  inverse  filtering.  A few 
restoration  examples  were  shown  to  illustrate  the  performance  of  the  new 
image  restoration  system. 
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APPENDIX 


In  the  Appendix,  we  provide  definitions  for  various  terms  used  in  this 
paper  such  as  discrete  space  Fourier  transform,  correlation,  power  spectrum 
and  energy  spectrum.  We  also  determine  "k",  the  normalization  constant  used 
in  equations  (4),  (6)  and  Figure  3. 


DEFINITIONS 


x(n^,n2):  a sample  sequence  of  a 2-D  stationary  random  field. 

f(nj,n«):  a finite  2-D  sequence  obtained  by  windowing  x(n^,n2)  with  a 
rectangular  window. 


f(n1,n2)  = 0 except  for  0 <_n1  ^N^-l  and  0 <_n2  ^N2~l 

N, -1  N„-l 


OO  00 


F (uij , o>2) : Z Z f(n1,n2)»e 


n^-00  n2=-°° 


-jw  n -jo) _n  1 2 

1 -e  1 Z-Z  £ f(n.,n,).e 

nj=0  n2=0 


-J“2n2 


discrete  space  Fourier  transform  of  ffnj,n2) 

2 

|F(u  ,u»  )|  : energy  spectrum  or  periodogram  of  f(n  ,n_) 

1 L A 4 

Rx^nl’n2^:  EtxCHj.fl^)  xCAj-nj,^-^)]  , correlation  of  x(nj,n2) 


OO  00 


P (w,  ,<0:  Z Z R (n, ,n0)*e 

x 1 2 xv  1 2' 

n =-<=°  n2=-°° 


-3u>1n1  -3u)2n2 


discrete  space  Fourier  transform  of  R (n,,n„)  or  power  spectrum 

x l i 


of  x(n1,n2). 
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Determination  of  "k" 


Given  f(n1>n2),  one  way  to  estimate  Rx(n1»n2^  is  by 


00  00 


VW  = N2N2  'i=_a  f(W  * ftVVW 


(Al) 


In  the  frequency  domain,  equation  (Al)  is  equivalent  to 


(A2) 


where  k = Nj-N2 


If  f(n1,n2)  is  obtained  by  a more  general  window  function  w(n1>n2),  a 
straight-forward  generalization  of  equations  (Al)  and  (A2)  are  given  by 


00  00 


R*<Vn2>  - - \S_  . * . 

*1  z 


(A3) 


A 2 

and  Px(w1»w2)  = £ • | F(o)1,o;2)  | 


(A4) 


00  00 


where  k = E E w (£,  ,£_) 
*!=-•  £2— • 1 


(A5) 
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